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ABSTRACT 

The Tully-Fisher relation of spiral galaxies shows notable dependence on morphological types, with 
earlier type spirals having systematically lower luminosity at fixed maximum rotation velocity V max . 
This decrement of luminosity is more significant in shorter wavelengths. By modeling the rotation 
curve and stellar population of different morphological type spiral galaxies in combination, we find 
the Vmax of spiral galaxies is weakly dependent on the morphological type, whereas the difference 
of the stellar population originating from the bulge disk composition effect mainly account for the 
morphological type dependence of the Tully-Fisher relation. 

Subject headings: galaxies: spiral - galaxies: stellar content - galaxies: kinematics and dynamics 



1. INTRODUCTION 

The Tully-Fisher relation (|Tullv fc Fisheri 119771 . here- 
after TFR) is an empirical relation between the absolute 
magnitude M and the maximum rotation velocity V msiyi 
of spiral galaxies, which is typically expressed as 



M = a Log 



V m 



200 kms 



0, 



(1) 



where the values of slope a and zero-point (3 are depen- 
dent on the photometric band of the M being measured. 

It has long been known that the TFR has a morpho- 
logical type dependence, with earlier type spiral s having 
systematically lower lum inosity at fixed V max (|Robertsl 
119781 iRubin et al.l fl980h . This luminosity difference is 
also waveband dependent, with larger offsets in shorter 
wavele ngths. For example, in / band, iGiovancll i et al.l 
(| 19971) found a 0.32 mag lower zero-point of the TFR 
for Sa/Sab galaxies and 0.10 mag lower for Sb galaxies 
than t he Sbc and later type spirals. In B band. TRusselll 
(|2004f ) found that Sb galaxies have a zero-point of 0.57 
mag l ower than Sc galaxies. In a recent study, IRussell 
(2008) compared the TFR-derived distances of nearby 
groups and clusters and found a mean difference of 0.19 
mag in H band between the Sb and Sc spiral galax- 
ies. Moreover, usin g a la rge sample of spiral galaxies, 
IMasters et all (|2006l l2008D found that the morphological 
dependence of the TFR is not only a shift of the zero- 
point, but even dependent on the luminosity in the way 
that the differences are more pronounced for more lumi- 
nous galaxies. 

The morphological dependence of the TFR originates 
from either the differences of the stellar population(M) 
or the disk dynamics (Vm ax ) or both. 

On one hand, it is well kn own that the colors ofear- 
lier type spirals are redder. iDevereux fe Youriel (|1991[ ) 
interpreted this color difference as originating from the 
bulge disk composition effect. The stellar population of 
bulges are typically older and more metal rich than that 
of the disks. A larger fraction of the bulge component in 



earlier type spirals naturally results in an redder color on 
average for the whol e galaxy. Howeve r with a bulge disk 
decomposed sample. [Kennicutt et al.l (|1994f ) studied and 
compared the star formation histories of only the disk 
component of different type spirals and found that the 
stellar population of the disks of later type spiral is also 
on average younger. 

On the other hand, although the rotation curves of 
spiral galaxies are prop osed to follow a universal shape 
(jPersic fc Saluccilll99ll iPersic et al.lll996ft . there is evi- 
dence showing that the rotation curves of early type spi- 
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massive Sa galaxies lie better in the well-defined TFR 
when using the asymptotic rotation velocity V^ymp in- 
stead of Vmax, implying a dependence of V max on the 
morphological type. 

On the theoretical side, the zero-point, slope and 
scatter of the observed TFR can all be well accommo- 
dated by the current disk formation model in the frame- 
work of the cold dark matter hierarchical cosmogonies 
dDalcanton et al.|[l997HMo et al.|[l998HMo fe MaofeoOOl ; 
iPizagno et al.ll2005h . However, the morphological depen- 
dence of the TFRs has not been probed in these studies 
because the bulge component in these models is typically 
neglected. Another limitation of these models is that the 
stellar populations have not been tackled with a physi- 
cal prescription, but with pre-determined mass-to-light 
ratios. 

In this study, we aimed to model the dynamics and 
stellar population of different type spiral galaxies in com- 
bination and try to find out which factor is the main con- 
tributor to the morphological dependence of the TFR. In 
specific, we will follow the dynamical model of Mo, Mao 
& White (1998, hereafter MMW) and extend it to in- 
clude a bulge component. For the stellar population, we 
will parameterize the star formation histories of the disks 
and bulges separately and then derive their mass-to-light 
ratios in dif ferent bands using the ste llar population syn- 
thesis code (|Bruzual fc Charlotj|2003[ ). 
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This paper is organized as follows. In Section 2, we 
describe our dynamical model of the spiral galaxies, in- 
cluding the disk, bulge and dark halo components re- 
spectively. In Section 3, we study the stellar properties 
of the disks and bulges with parameterized star forma- 
tion histories. We compare our model predictions with 
the observational data in Section 4. We make discussions 
on the uncertainties of our results in Section 5 and finally 
give a brief summary in Section 6. 

2. DYNAMICAL MODEL OF SPIRAL GALAXIES 

Our modeling of the dynamics of spiral galaxies fol- 
lows and simplifies the disk formation model of MMW, 
but with more attention on the bulge contribution. The 
interested reader is referred to MMW for detail. Here we 
repeat the essentials related to our study. In this model, 
initially, the gas and dark matter are uniformly mixed 
in a virialized halo. As a result of dissipative and ra- 
diative cooling, the gas gradually cools down and settles 
into a disc structure due to the conservation of angular 
momentum. 

We define that the stellar mass of the galaxy finally 
formed is M* and the fraction of this mass to the initial 
halo mass Mh is m s = M^/Mh- We express the bulge 
fraction of the formed galaxy being /t>, so that the masses 
of the bulge and disk are Mb — and M d = (1 — 

/b)M* respectively. 

The rotation curve of a spiral galaxy is contributed by 
three dynamical components: halo, bulge and disk. We 
show and discuss our assumptions on each term below. 

2.1. halo 

The iV-body simulations show that the collapsed and 
viri alized dark halos follow a univers al density NFW pro- 
file (jNavarro. Frenk. fc Whitelll996t ). 



p(r) 



V h 2 



r/r s 



4irGr 2 [ln(l + c) - c/(l + c)] (r/r s + l) 2 



(2) 



where r s is a scale radius, G is the gravitational constant, 
Vh is the circular velocity and c is the concentration pa- 
rameter. The concentration c is defined as c = r2oo/f s , 
where r2oo is the virial radius 1 of the halo. For a given 
cosmology, the viral radius r2oo, halo mass Mh and cir- 
cular velocity Vh are related by 



Vh 



^200 



ioH(zy 



M h 



v h 3 



IQGH(z) 



(3) 



where H(z) is the Hubble constant at redshift z. 

In this study, we use the concordance ACDM cos- 
mology model, with H = 70 kms~ 1 Kpc~ 1 , £lo = 0.3, 
Qa = 0.7 and baryon mass density fig = 0.04. 

With the assembling of the baryons into disk and bulge, 
the gravitational effect from the bulge and disk changes 
the initial halo mass distribution through contraction. 
We follow MMW and use the adiabatic contraction as- 
sumption to analysis this effect (see MMW for detail). 

For a halo with given circular velocity Vh , the concen- 
tration c is the only parameter to be quantified. At given 
redshift z, the concentration param eter c is mainly cor 
related with the mass of the halos jNavarro et al 



lycor- 
] tl99ll 



1 Here, r2oo is the radius within which the mean mass density is 
200 times the cosmic critical density p cr it- 



iBullock et alj 120011 ) . We adopt a simple parametriza- 
tion of the conce ntration c at redshift zero as that in 
IShen et all ((20(51) : 



c = 8.5 



V h 



100 kms- 



2.2. bulge 



-1/3 



(4) 



We assume that the bulge component of spiral galax- 
ies has a spherical mass distribution and the mass den- 
sity profile pb{r) follows a Hernquist profile (Hernquist 
1990), whose projection approximates the classical i? 1 / 4 
surface brightness profile of elliptical galaxies. The Pb(r) 
is expressed as 



, v M b a 



2-7T r (r + a) 3 



(5) 



where Mb is the total bulge mass, a is the bulge scale 
radius, which is correlated with the eff ect ive (half- light ) 
radius R e in the way R e « 1.82a[Equ. 38 of iHernauisti 
(1990)]. To establish the dynamics of a bulge with mass 
Mb, we need to know the scale radius a or effective radius 
R e . We assume that the bulges follow the ob served size- 
mass i (R — M) relation of elliptical galaxies (IS hen et al.l 
2003). 



E e (Kpc) = 3.47 x 10" 



'M h 



0.56 



(6) 



2.3. disk 



The surface brightness profile of spiral disks typically 
follow an exponential profile. Here, we assume the disk 
surface mass density profile is also exponential, 



Pd(r) = p exp(-r/R d ) 



(7) 



where po is the central surface mass density and Rd is 
the scale-length, po and Rd are related to the total mass 
of the disk Md through Md = 2ttpqR^. 

Since the observed size of the disk in a spiral galaxy is 
affected by its central bulge, we do not take the observed 
R — M relation of spiral galaxies as that in the model. 
Following MMW, we assume that the disk size Rd is de- 
termined by its initial angular momentum J, which can 
be parameterized by a spin parameter A 



A = J\E\ 1/2 G- l M~ 5/2 , 



(8) 



where E is the total energy of the halo. With reason- 
able assumptions that the specific angular momentum 
per particle of the baryons is the same as the dark mat- 
ter and there is no angular momentum transfer among 
different components, the angular momentum of the disk 
finally formed is therefore 3d = mdJ- Here, md is the 
fraction of baryons settling into the disk and is equal to 
m s (l — fb). The scale-length of the disk Rd then equals 



Rd 



±\r WQ f c 



-0.5 



f R (X,c,m d ) 



(9) 



where f c is a factor coming from the halo density profile 
and only dependent on concentration c (see Equ. 23 of 
MMW) , fn is a factor coming from both the halo density 
profile and the gravitational effects of disk and bulge(see 
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Equ. 29 of MMW). It is shown that the size of spiral 
galaxies predicte d from A is in excellent agreement with 
the observations (|Sver et al.lll999t iShen et al.ll2003f) . 

Numerical simulations have found that the spin param- 
eter A of the dark matter halo follows a log-normal distri- 
bution with median A ~ 0.04 and scatter o~\ n \ ~ 0.4 and 
this distribution is quite independen t of the cosmology 
and halo properties (jBett et al.l 12007). For the purpose 
of simplicity, we do not consider the scatter of the an- 
gular momentum distribution and take the median value 
A = 0.04 for all disks(see more discussions in Section l5TTj) . 

2.4. rotation curve 

With the mass distribution of the bulge, halo and disk 
components determined through the model parameters 
shown above, the rotation curve of a spiral galaxy can 
be determined by the sum in quadrature of contributions 
from these three terms: 



V 2 (r 



= c M DM (r) + M b (r) + y2 



Vj(r). (10) 



During the calculation of the term V d 2 contributed from 
the disk component , the flattened geometry n eeds to be 
taken into account (jBinnev &: Trem ainc 1983). 

To fully determine the rotation curve of a given spiral 
galaxy with stellar mass M*, we still need to know the 
relative proportions of these three components (bulge, 
disk and dark halo). This could be done by settling the 
two other parameters, i.e. the baryon fraction m s and 
the bulge-disk-ratio B/D 2 . We set the baryon fraction 
m s to be a f unction of the hal o mass (in terms of circular 
velocity V>.) (|Shen et al.ll2003h . 



0.13 



1 + (Vfc/150 kmt 



(11) 



This assumption is based on the consideration that 
gas outflow process in low mass halo is efficient, 
due to their shallow potential well, and a large 
fraction of baryons is blown out of the dark halo, 
while the baryon fraction of high mass halos gradu- 
ally approach the cosmic baryon fraction £Ib/Qq ~ 
0.13 (jWhite fc Frenkl fl99lf) . The bulge-disk-ratio B/D 
is taken to be 0.5/0.3/0.1 for model Sa/Sb /Sc galaxies 
respectivelv ljSimien fc de Vaucouleurslll986l ). 

We show three examples of resulting rotation curves 
for three morphological types(Sa/Sb/Sc) in Fig. [TJ All 
three model galaxies are chosen to have the same circular 
velocity Vh = 120 km s -1 and thus the same stellar mass 
M = 3 x 1O 1O M . For the earlier type spiral, as we can 
see from the figure, due to the larger bulge component, 
our model predicts a steeper increase of rotation velocity 
in the inner region, which brings on the appearance of 
Knax at smaller radius, and a systematically larger V max . 

However, the value of V max only weakly correlates with 
the model type. The V max of the Sa(Sb) spiral is only 
1.5(1.0) percent higher than the corresponding Sc spiral. 
This is because the relatively smaller bulge fraction of 
the later type spiral is always compensated by its larger 
disk fraction (see the components of rotation curves in 
Fig. ©. 

2 B/D = f b /(l - f b ) 
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Fig. 1. — The rotation curves of the model galaxies with Vh = 
120 kms -1 . The solid lines show the rotation curve while the dot- 
dashed, dotted and dashed lines show the contributions from disk, 
halo and bulge components respectively. The green, red and blue 
lines represent the model Sa, Sb and Sc galaxies respectively. 

3. STELLAR POPULATIONS OF SPIRAL GALAXIES 

Besides the different dynamics, the bulge and disk also 
show different stellar populations. The bulges are typi- 
cally old and show little recent star formation, while the 
star formation time-scales of the disks are long and there 
is still ongoing star formation throughout the disk at the 
present day. 

In this study, we parameter the stellar population of 
the bulge components with a single stellar population 
with age of 10G year(Gy). For the disks, following the 
usual convention, we parameterize their star formation 
histories(SFH) with an exponential function 



SFR(t) cx exp(-i/r) 



(12) 



where r is the time scale of the SFH. This simple analyt- 
ical expression can parameterize different kinds of SFHs. 
For r approaching 0, Equ. (Ti"2"|) represents a single stellar 
population. When r tends to infinity, the star formation 
rate is a constant along the history. The age of the disks 
is also set to 10 Gy. 

Kennicutt et al. (1994, hereafter K94) parameterized 
the SFH of the spiral disks with a b parameter, which is 
defined as 

SFR , . 

where SFR is the star formation rate today and < 
SFR > is the average star formation rate in past. 
They estimated the b parameter from UBV colors 
and Ha fluxes for different morphological type spirals 
and found that the Sa/Sb/Sc disks have typical b w 
0.12/0.33/0.84 and variation of about 0.05/0.15/0.20 re- 
spectively. 3 With the SFHs parameterized by Equ. [T2l 

3 The b values are taken from Table 4 of K94, where the sub- 
types include Sa/Sab/Sb/Sbc/Sc/Scd spirals. We average the b 
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as a prior, these b values correspond to the r with me- 
dian 3/5/30 Gy and la range [2,4]/[4,7]/[10,oo] Gy for 
Sa/Sb/Sc spirals respectively. 4 

Based on the SFHs assumed for disk and bulge com- 
ponents, we use t he new version of th e stella r population 
synthesis code of iBruzual fc Charlotl (|2003f ). denoted as 
CB07, to calculate their mass-to-light ratios in different 
wavebands. The metalicity of the stellar population is 
set to be the solar val ue while t he ste llar initial mass 
function is taken from [Chabricr (|2003l ) with lower and 
upper mass limit being 0.1M Q and 100M© respectively. 

Finally, with the predetermined stellar mass m s and 
B/D for each model galaxy, we calculate its absolute 
magnitudes in different optical wavebands. 

4. MODEL PREDICTIONS 

In this section, we show our model predicted morpho- 
logically dependent TFRs and compare them with the 
observational results. For each morphological type, we 
build 11 model galaxies with halo circular velocity Vh 
in the range from 100 to 200 kins" 1 with an interval 
of 10 kms -1 . With the baryon fraction characterized 
by Equ. [11] the range of the stellar mass of our model 
galaxies is thus from 1.4 x 10 10 to 2.3 x 10 n M Q . 

In our model, there is only one free parameter, i.e. 
the disk star formation time scale r in Equ. [T2l We 
first set the r values of different type spirals to be these 
suggested by K94, i.e. the median r = 3/5/30 Gy and its 
la range being [2,4]/[4,7]/[10,oo] Gy for Sa/Sb/Sc spirals 
respectively. We refer to this model setting as K94 model 
below. 

4.1. I-band TFRs 

Lots of recent observational TFR studies a re in 
the near-infrare d J bandfe.g. IGiovanelli et al.l 119971 : 
iDale et al.l[l999l : iVogt et al.l I2004H . which takes the ad- 
vantage of better detector photometry than H band and 
less scatter in TFR than blue B band. More recently, 
iSpringob et alj (|2007h compiled and published a new gen- 
eration of a homogeneous galaxy catalog, referred to as 
the SFI++, which contains ~ 5000 spirals suitable for 
TFR study in / band. Based on the SFI++ catalog, the 
T FRs of Sa, Sb a n d Sc type spirals have been presented 
bv iMasters et al.l (|2006l hereafter M06). Here, we first 
show our model predicted / band TFRs for different type 
galaxies and make comparisons with the results of M06. 

The K94 model predicted / band TFRs of Sc, Sb and 
Sa spirals are shown in the top left, top right and bot- 
tom left panels of Fig. [2] respectively. The model galax- 
ies are represented by the triangles(Sc), squares(Sb) and 
stars(Sa) while the fitted TFRs are shown as the solid 
lines in each panel. The observed / band TFRs of dif- 
ferent type spirals of M06 are shown as the dashed lines 
for comparison. The dotted lines in each panel show the 
ranges of model predicted TFRs when the disk star for- 
mation time scale r varies in its la range. Inside its 

values of Sa and Sab sub-types as the Sa type. The b value of 
Sb type is kept. For Sc type, we average the b values of Sbc and 
Sc/Scd sub-types. The variations of b are estimated from Fig. 6 of 
K94. 

4 Here, the la range denotes the 32th to 68th percentiles of the 
r distribution and is calculated from the estimated variation of b. 
For Sc spirals, during the calculation of the range of r, we force 
the upper limit of b to be 1, corresponding to r ~ oo, which g ives 
the youngest average stellar population as that can do by Equ. 1121 
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Fig. 2. — The model predicted I band TFRs of different mor- 
phological type spiral galaxies. The triangles, squares and stars 
represent the model Sc, Sb and Sa galaxies in top left, top right 
and bottom left panels respectively. The solid lines show the best 
TFR fittings of model galaxies while the observed TFRs of M06 
are shown as the dashed lines for comparison. The dotted lines in 
each panel show the range of the model predicted TFRs when the 
star formation time scale r varies in its ltr range (as that denoted 
in the brackets). 

range, a larger r results a larger zero-point of the TFR, 
i.e. galaxies will be brighter at given Vmax for larger t. 

Sc type spirals (include Sbc/Sc/Scd) are the domi- 
nant morphological type in TFR studies 5 . Lots of 
the TFR studies use Sc spirals as the reference and 
make morpholog i cal ty p e corrections on o t her types, e.g. 
IGiovanelli et all (|1997l ): IDale et all (|1999f ): IZiegler et al.l 
(2002). For Sc spirals, as we can see, by setting the model 
parameters to be either the typical values from numeri- 
cal simulations(e.g. A and c) or being constrained from 
other observations(e.g. r ~ 30 Gy) without any further 
fine-tuning, our model predicted / band TFR is in excel- 
lent agreement with the observations. 

However, with default settings of model parameters, 
our model predicted TFRs of Sb and Sa spirals are not 
well consistent with the observations of M06. The pre- 
dicted TFR of Sb spirals(r=5 Gy) have a systematically 
larger zero-point than observations, i.e. the model pre- 
dicted luminosity is too bright. When r reaches its upper 
range 7 Gy, the discrepancy becomes smaller, but still 
be far from consistent with observations. We will discuss 
this discrepancy in more detail in next section. 

For Sa spirals, the slope of the observed TFR is signif- 
icantly shallower t han our model predicti ons. However, 
as that discussed in IGiovanelli et all (|1997t ) and M06, the 
disagreement of the TFR slope of the Sa spirals should 
be interpreted in caution because the incompleteness of 
the Sa spirals is large, especially in the low luminos- 
ity end, which will artificially bias the slope to lower 
value. On the other hand, our modeling on the B/D of 

5 The later type spirals are more likely to be strong HI or Ha 
emitters so that the measurement of the rotation width is easier. 
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Sa spirals is a constant 0.5, which might be too simpli- 
fied. In the study of M06, all the SO/Sa/Sab galaxies 
are grouped as Sa type. Among three sub-types, the ear- 
lier type galaxies, e.g. SO, having systematically larger 
B/D, are also systematically brighter. Therefore, from 
low V max to high Vmax galaxies, there is actually a sys- 
tematical change of B/D, which will also shallow the 
slope of TFR. By introducing a systematical change of 
B/D as function of Vmax, we can in principle reproduce 
the slope of Sa spirals as observed. However, we prefer 
not to tune our model to reproduce this relation with a 
cost of introducing more uncertainties. 

Comparing with the uncertainties in the modeling of 
the TFR slope, the zero-point is quite robust. More- 
over, many studies of the morphological dependence of 
the TFRs only repor t a glo b al shift of the zero- point, 
e.g. IGiovanelli etaTl (Il997h : IRusselll IpOOi l200l (see 
next section). Therefore, in the following section, we pa- 
rameterize the morphological dependence of the TFRs 
only with the shift of the zero-point AM. 

4.2. morphology dependent TFRs in different bands 

We choose the zero-point shift AM of the TFRs at 
Vmax = 200 km s -1 , where the corresponding absolute 
magnitude is roughly the M» of the luminosity function 
of spiral galaxies and where the number counts of galax- 
ies in a flux-limited sample normally peaks. In this case, 
AM = A/3 of Equ. 1. 

We show the AM between Sb and Sc spirals in the 
top panel of Fig. [3l whereas the AM between Sa and 
Sc spirals is shown in the bottom panel. The model pre- 
dicted AM in different wavebands are connected as a 
function of their effective wavelengths with lines. The 
wavebands include U,B,V,R,I,J,H,K and the wave- 
length ranges from 0.33 to 2.2 micron. The observed dif- 
ference of AM at Vmax = 200 km s -1 in different bands 
found in the literature are plotted against their effec- 
tive wavelengths. Results from different authors are la- 
beled with different symbol types. The references and 
details of these morphology dependent TFRs are listed 
in Table 1. A few other studies in the l iterature are 
not quoted in Table 1 and Fig El including iRubin et all 
(fl985h. IGiovanelli et alJ (Il997h and the TFR of Sb spi- 
rals of ISandaed (|2000f) T IRubin et all (fl985l) found the 
magnitude differences AM between Sa and Sc spirals 
are as high as 2 mag in B band and 1 mag in H band, 
whose number of each type galaxies is limited(< 20) and 
without error estimation q uoted for the paramete rs of 
TFRs. The / band data of IGiovanelli et all (|1997f) has 
been expanded and re-analyzed by M06 and their re- 
sults are genera lly consistent. The TFR of Sb spirals in 
Sandage (2000) is arbitrarily neglected, where its zero- 
point is even lower than Sc spirals, i.e. AAf(Sb-Sc)< 0, 
contradictory to most of the other studies. 

The solid lines in Fig. [3] show predictions from the 
K94 model, i.e. the model with r = 3, 5,30Gy for 
Sa,Sb,Sc spirals respectively. The shadowed regions show 
the ranges of model predicted AM when r varies in its 
lcr range. As we can see, the K94 model over-predicts the 
differences of the zero-points significantly for both Sb/Sc 
and Sa/Sc pairs. The average observed zero-point differ- 
ences in near-infrared bands are AM(Sb-Sc)~0.15 mag 
and AM(Sa-Sc)~0.20mag, whereas the model predic- 
tion is as high as AM(Sb-Sc)~0.45 mag and AM(Sa- 



Sc)~0.60mag. That means, if the differences of the stel- 
lar populations of the disks of different type spirals are 
as that suggested by K94, the shift of the zero-point of 
the TFRs would be much larger than observed. 

With the disk star formation time scale r approach- 
ing its upper range (the lower boundary of the shadowed 
region), the discrepancy between the model predicted 
AM and the observations becomes smaller. This result 
prompts us to consider an extreme case, r = 30Gy for 
both Sa and Sb disks. In this case, the stellar pop- 
ulations of the disks of different type spirals are the 
same, while the different global colors of different type 
spirals only originate from the bulge disk composition 
effect. This scena ri o is consistent with the results of 
IDevereux fc Youngl (|199lD . To distinguish it from the 
K94 model, we refer to this scenario(r = 30 Gy for all 
type disks) as the 'composition' model below. 

For the 'composition' model, the model predicted AM 
as function of the effective wavelength is shown as the 
dotted lines in Fig. [3] Surprisingly, the predicted AM 
from this model, especially in near-infrared bands, is gen- 
erally consistent with observations, cither for Sb/Sc or 
Sa/Sc pairs. The predicted AM in B band is smaller 
than that observed. A possible solution to this discrep- 
ancy is the dust extinction which has not been taken into 
account in our model. If the face-on dust extinction of 
earlier type spirals is more significant, the AM in blue 
band will be larger than now we predicted(see detailed 
discussions in Section l5~3]) . 

Finally, it is worth mentioning that the AM between 
different spiral types is not contributed by the stellar 
population alone but composed by both the dynamics 
and stellar population. As we have shown in section [2T4l 
for the galaxies with the same stellar mass 3 x 1O 1O M0, 
the Vmax of a Sa(Sb) spiral is about 1.5(1.0) percent 
higher than Sc spiral, which corresponds to the AM 
~ 0.04(0.03) mag, independent of the wavebands. This 
contribution is shown as two dashed horizontal lines in 
the top and bottom panels of Fig. [3] However, as we can 
see, the stellar population is still the dominant contrib- 
utor to the morphological dependence of the TFR. 

5. DISCUSSION 

There are some uncertainties in our modeling of both 
dynamics and stellar population of the spiral galaxies, 
e.g. the scatters of model parameter, the diverse bulge 
properties and the internal dust extinction etc. We dis- 
cuss these issues below. 

5.1. scatters of model parameter 

In our modeling, we only take the typical values for 
each model parameter and do not consider their scatters 
except the key parameter r. However, all these model 
parameters show scatters as sug gested by either numeric 
simula t ions or obse r vation s, e.g. IMo fc Mad (|2000[ ); [jmgl 
(120001) ; iShen et alJ (|2002D . If the scatters of model pa- 
rameters are independent variables and do not correlate 
with the morphological types(as suggested in our model), 
these scatters will only make contributions to the scat- 
ters of the p redicted TFRs and will not affect a ny of our 
conclusions (|Mo fc Maoll2000l : IShen et al.|[2002l ). 

However, the morphological type may correlate with 
the spin value of a galaxy. A smaller A may trigger a 
formation of larger bulge easier through disk instability 
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Fig. 3. — The morphology dependence of the TFR in different 
wavebands. The solid and dotted curves show the predicted shift 
of the zero-point of the TFRs, AM, as function of wavelength from 
the K94 model and 'composition' model respectively (see text). The 
dashed horizontal lines show the contributions to AM from the 
Vmax(see text). The upper panel shows AM between Sb and Sc 
spirals, while the AM between Sa and Sc is shown in the bottom 
panel. The observed AM in B, I, J, H and K band are labeled for 
comparison. 

(|Shen et al.ll2003l ), so that induce an earlier type spiral. 
In other words, the earlier type spirals would be biased 
to the galaxies with systematically lower A when we con- 
sidered the scatter of A. In this case, the resulted V max 
of earli er type spirals wil l be even larger because of the 
lower A ([Mo fc Maoll2000D and thus the room left for the 
stellar population in AM is less. 

5.2. pseudo and classical bulges 

In this study, we have treated all the bulges as scaled 
elliptical, i.e. the classical bugles, for simplicity. 

The current view on the bulges is quite complex. There 
are mainly two types of bulge, i.e. the classical bulge and 
the pseudo-bulge, which may be originated from differ- 
ent physical formation processes ([Fisher fc Drorvl l2Q08t 
iFisher et al.ll2009t iGadottil [20091) . The classical bulges 
are similar to scaled elliptical galaxies, with their surface 
brightness following classical i? 1 / 4 profile and colors be- 
ing red. The pseudo-bulges tend to show younger stellar 
population and their density profiles can be character- 
ized by a Sersic profile with Sersic index n significantly 
smaller than 4. The classical bulges typically appear in 
earlier type spirals and their masses are high. The major- 
ity of th e stellar mass o f the bulges is in classical bulges 
[~ 90%. [Gadottil ((20091) ]. 

If the bulges of late type(e.g. Sc) spirals are pseudo- 
bulges, the difference of stellar population between early 
and late type spirals will be even larger because the 
pseudo-bulges are bluer than classical bulges. This will 
increase the predicted AM between the early and late 
type spirals and decrease the allowed differences between 
their disks again. 



5.3. dust extinction 

In our modeling, we have not considered the effect of 
internal dust extinction, which might be different for dif- 
ferent type spirals. 

Fortunately, the internal dust extinction has at least 
been partly corr ected in most of th e observational 
TFR s t udies , e.g. IGiovanelli et all (|1997t ): IMasters et all 
(|2006l 120081 ). In these studies, they all made the 
corrections of the internal dust extinct ion using the 
parametrization of IGiovanelli et aTl ([1994D . 

A = 1 Log(a/b) (14) 

where 7 is a waveband dependent coefficient and a/b is 
the observed axis ratio indicating the inclination angle 
of the spiral disk. This parametrization of the dust ex- 
tinction has accounted the extinction from the geometry 
effect, i.e. it has corrected all the spiral disk to face- 
on viewing and assumed zero dust extinction for face-on 
case(A ~ for a/b ~ 1). 

However, the internal d ust extinction of face-on spi- 
rals m ay not be negligible (|Shao et al.ll2007t iDriver et al.1 
120071) . The inner parts of the spiral disks are very 
lik ely to be optical thic k even when viewed in face- 
on (iGiovanenriESIIH^). The dependence of this face- 
on extinction on the morphological type is complicate. 
On one hand, the later type spiral disks are gas riche r 
and so that may contain more dust (|Stevens et al.| [2005) . 
On the other hand, the earlier type spirals have larger 
bulge component and this component suffers more from 
the dust extinctio n because the dust layer is peaked i n 
the central region ([Tuffs et al J 12004 ; IDriver et alll2007l) . 
Counteracting these two effects, the global dust extinc- 
tion may only be weakly depende nt on the morpholog- 
ical ty pes. A very recent study of IMunoz-Mateos et al.l 
(2009) based on the multi- wavelength data shows that 
there is a weak global trend of the internal dust extinc- 
tion, Sb>Sa>Sc. 

If the face-on dust extinction of early type spirals 
were more significant considering their more concen- 
trated light, the smaller AM predicted by the 'compo- 
sition' model in B band could be explained since the 
extinction in B band is 5 times more larger than K 
band given the typical extin ction curve of normal spi- 
ral galaxies( |Shao et all 120071 ). Moreover, if we took the 
higher dust extinction for the early type spirals into our 
model, the model predicted AM would be larger, thus 
the stellar populations of different type disks even could 
not be very different as that suggested by K94. 

6. SUMMARY 

In this paper, we have studied the morphological 
dependent TFRs by modeling the dynamics and stel- 
lar population of spiral galaxies in combination. We 
model the dynamics of spiral galaxies by following 
the disk formation model of MMW and extending it 
to include a bulge component. We parameterize the 
SFHs of the bulge and disk separately and derive their 
mass-to-light ratios with population synthesis code of 
IBruzual fc Chariot] ([2003D . Our model reproduces the 
observed / band TFR of the Sc spirals very well without 
any free parameters. Our model shows that the morpho- 
logical dependence of the TFRs is mainly contributed 
by the stellar population through bulge disk composi- 
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tion effect, although the effect from the dynamics is not 
negligible. 

The effects of the dynamics and stellar population on 
the morphological dependence of the TFR are different. 
The shift of the Vmax (or luminosity) , caused by dynam- 
ics, shows a characteristic that is independent of the 
wavebands. On the other hand, the shift of the zero- 
point of TFR is a function of the wavebands, when the 
morphological dependence is caused by the stellar popu- 
lation. Therefore, the morphological dependence of the 
TFR in different bands is an effective way to constrain 
the dynamics and stellar populations of spiral galaxies. 
Our results show that the observed morphological depen- 



dent TFRs are consistent with a scenario that the stellar 
population of different type spirals are only different in 
the bulge disk composition on average. 
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TABLE 1 The observed morphological dependent TFRs in literatures. All 
the TFRs are re-shaped in the form of Equ. (1). AM is the shift of the 
zero-point /3 relative to the Sc spirals in the same study. 
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